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Abstract 

We calculate the damping of excitations due to four-fermionic interaction 
in the case of two-dimensional superconductor with nodes in the spectrum. 
At zero temperature and low frequencies it reveals gapless lo 3 behavior at 
the nodal points. With the frequency increasing the crossover to the normal- 
state regimes appears. At high frequencies the damping strongly depends 
on details of a normal-state spectrum parametrization. Two important par- 
ticular cases such as the models of almost free and tight-binding electrons 
are studied explicitly and the characteristic scales are expressed through the 
model-free parameters of the spectrum at the nodal points. The possibility 
of crossover in temperature dependence of damping in the superconducting 
phase is discussed. 
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I. INTRODUCTION 



An assumption that the pairing state in the high-temperature superconductors is the 
ci-wave-like one is widely discussed nowadays. Q Although up to now it apparently cannot 
be treated as the uniquely established fact there exist a large body of experimental data 
which supports it. Such a pairing would result in presence of "nodal" points in the su- 
perconducting excitations spectrum in whose vicinities it has gapless phonon-like character. 
It should essentially modify the damping behavior comparing to the s-wave pairing state. 

The quasiparticle damping owing to impurities has been extensively studied in literature. 
U However, this is not the only channel, though, maybe, the more interesting one. The 
excitations scattering by phonons and also by each others always contribute to the damping, 
as well. 

In the present paper we evaluate the quasiparticle damping in two-dimensional (2D) 
superconductors with d-wave-like pairing appeared due to their four-fermionic interaction. 
It may be the on-site Coulomb repulsion in the small-f/ Hubbard model, or the residual part 
of the interaction which provides the pairing in the strong correlations scheme. 

In the BCS case, the damping at zero temperature has the threshold at u = 3A. We 
show that for the gap function with nodes it demonstrates gapless u 3 low-frequency behavior 
at the nodal points. The results obtained for high frequencies depend on details of Fermi 
surface (FS) parametrization. The damping for the circular FS will increase quadratically 
while for the spectrum with Van-Hove singularities the additional linear regime appears. 
Being generalized onto finite temperatures, this results allow us to conclude about existence 
of crossovers in temperature dependence of damping in the superconducting state. 

The outline of this paper is as follows. In Sect. II we obtain general expression for 
damping of superconducting excitations due to their interaction. In Sect. Ill we evaluate 
it for the case of 2D superconductor with nodes in the spectrum. Since the high-frequency 
behavior of damping depends strongly on FS parametrization, in the consequent subsections 
we analyze three principal cases: (i) the case of circular FS and also the case of FS of varying 
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curvature when the gap node intersects it (ii) at the almost flat piece or (iii) in the vicinity 
of saddle-point of the spectrum. At last, in Sect. IV we generalize and discuss the results 
obtained. 



II. GENERAL FORMALISM 



It is convenient to represent the Hamiltonian Ji = H + 'Hi in the Nambu notation 

n o = J2 *k (fe + A k r- + A*r+)M/ k , (1) 

k 

Wi = \ £ Vi c+q)k ,_ q;kik ,(*t +q T3* k )(*t_ q r 3 * k ,). 

Z k,k',q 

Here Vk 1 k 2 ;k 3 k 4 is the interaction, = (c k j , c_ k jj stands for the spinor with c kcr being 
the fermionic creation operator, £ k and A k are the normal-state conduction band and the 
gap function (GF), respectively, and = \ (ti ± ir 2 ), where denote the Pauli matrices. 



Let us define the bare matrix Green function as follows 

/ 



£(0) 



C(o) _p(o) \ 
F (o)+ G (o)+ ) 

where normal G and anomalous F Green functions are written as 



V 



(2) 



^(0) / • \ i^n + tik 



F (0) ( 



(3) 



Here iu n = 2iri(n + |)T is the imaginary Matsubara frequency and the excitation spectrum 
is given by e k = [£ 2 + | A k | 2 ] 1/2 - 

Within this notation the Dyson equation obtains the form G^ 1 = G^^ 1 — E, with the 
diagonal self-energy operator E = diag (E, E + ). 

The first nonvanishing contribution to E is given by 



E(k,?W n ) = T 2 ^l,2;3,oV3,0;l,2 5l,2;3,0 



1,2,3 



X 



(4) 
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where the indices 1, 2, 3 label the intermediate frequencies and momenta and corresponds 
to the external one. 

After the analytic continuation iu n — > uj + i5 for the imaginary part of the self-energy at 
T = and uo > we get 

UJ LO — X 

ImS(k,^) = -7T^ Wlfe^q+p^^p^p^q+p^ 

wo 

x 5 (x-Eq) 5 (y-Sp) 5 (a;+y-u;+£ q+ p_ k ) 

w r 2 2 2 * * 2 i /c-\ 

x [vVq+P-k-^q^VVVW-kh ( 5 ) 
where the coefficients of Bogolyubov transform are defined by 

III. CALCULATIONS OF DAMPING 

In this section we will evaluate Eq. (5) in the case of 2D superconductor with unconven- 
tional pairing assuming different model parametrizations for the Fermi surface. 

We shall suppose that the gap function has nodal lines | A q | = in the 2D Brillouin zone. 
Let these lines intersect Fermi surface at a set of the points q = q n . Then in the vicinities 
of such points one can write 

£ q ^ v F q, A q ~ v 9 q; 

(7) 



<9q 



<9A q 

V 9 = 



q=0 



9q 



q=0 



where q = q — q„ and the index n labels the nodal point. We see that the excitation 
spectrum near these points has gapless linear character. 

In general, the quasiparticle damping due to their interaction has to reveal threshold 
feature in frequency in the superconducting state. The threshold corresponds to the minimal 
energy of the excitation needed to decay into three quasiparticles in the intermediate state. 



For the gap function with nodes however it disappears at the nodal points of the spectrum. 
Consequently, the damping at this points appear to be gapless. Taking it into account, we 
restrict ourselves in treatment of damping in the vicinities of nodal points only. 
The straightforward calculations at u — > yield 

ImE(k n ,w) ~ -X n ^ 

n / j I / 7 11 h h, 

Vn;lk = . li, m ^k,q+p-k;q,p^q,p;k,q+p-k, 

where the indices /, k run over all the nodal points. Here Tik are some dimensionless positive 
factors. The diagonal terms of the matrix T (i.e., referring to a single nodal point) are always 
numbers of order of unity while the off-diagonal ones could have additional smallness for 
particular FS parametrizations via their dependence on spectrum parameters. The principal 
contribution to the damping at low frequencies comes from excitations with momenta lying 
in the small elliptic vicinities of nodal points (see Fig. |1|) wherein linear approximation for 
the spectrum works well. 

Note that we obtain the universal cubic frequency dependence of the excitations damping 
at the nodal points without appealing to any fixed parametrization for spectrum. By the 
way, such dependence means the stability of nodes with respect to this kind of interaction. 

At the same time, some features of the model may, in principle, lead to different power 
laws. In particular, the damping would increase with frequency slower than a; 3 if at least one 
of spectrum zeroes (i) appeared due to tangency of FS and GF node, or (ii) coincided with 
van- Hove singularity in the band structure and/or (iii) with cusp- or discontinuity-point of 
the nodal line of the gap function. All these cases are obviously quite exotic. 

When the frequency increases, the integration in Eq. (5) cannot be restricted to the 
nodal vicinities and such characteristics as general shape of FS, curvature of FS at the 
nodal points, existence of van-Hove singularities in the spectrum, etc., should be taken into 
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account. One can distinguish three qualitatively different cases: (i) the circular FS and 
also the FS of varying curvature with the nodal points located (ii) far from or (iii) at the 
sharp curvature maximum. We shall illustrate the above cases by examples which would 
allow us to examine in detail the intermediate asymptotics and to express the corresponding 
crossover scales through the parameters of the models used. 



A. Circular Fermi surface 



Putting the lattice parameter to be unity, we write the circular FS as follows: £ q = 
q 2 /2m — Sp. For definiteness, we choose the d^.^-wave gap function, A q = A/2 [cosg x — 
cosq y ]. HQ In line with Eq. (8) we get: 



/ V \ uj 3 A 2 
ImE(k n ,^)~- — — , , (9) 

\e F J A z e F 

where V 2 ~ msixV 2 . ls . The principal contribution to ImS is obtained when all the momenta 

p, q and k n are almost parallel and belong to the vicinities of the nodal points (see Fig. [I]) . 

At higher frequencies A 2 /ef «C uj <C e F we find 

/ y n 2 2 
lmZ(k n ,u)--[—) —In 



UJ Ef £f\ 

U)JUi A^^J 



(10) 



^£p J £p 

To obtain Eq. (10), the curvature of FS was taken into account (see Appendix). The 
crossover from the cubic dependence to the quadratic one appears at uj ~ A 2 /ep. It can 
be qualitatively explained on the following way. One can show that, when uj A 2 /sf, 
the areas of integration defined by arguments of 5-functions in Eq. (5) are located near the 
nodal points (see Fig. [I]) where the spectrum can be well approximated by ellipses. When 
uj ~ A 2 /ef, the areas of integration propagate along the whole Fermi surface so that not only 
the vicinities of the nodal points principally contribute to ImS. For A 2 jep uj < A, we 
obtain the quadratic intermediate asymptotic regime with unconventional logarithmic factor 
of the form Iii(uj£f/A 2 ) which stems from integration over the FS. At frequencies higher 
then the gap value A < w < e f information about the superconducting state disappears, 



6 



and the conventional 2D Fermi-liquid result JF| is restored. Note that principal contributions 
to ImE in the latter case again come from almost parallel k, p and q. 



B. Tight-binding model, d x 2_ y 2-wave pairing 



Now we will demonstrate how the deformation of FS affects the damping. As before 
we shall consider d x 2_ y 2 gap function. Let us adopt the simplest tight-banding form for 
the normal-state dispersion: £ q = — It [ cos q x + cos q y ] + fi. Here t is the nearest-neighbors 
hopping integral and fi is the chemical potential measured from the half-filling. Note that for 
small enough fi -C t Fermi surface has long almost flat pieces near the points of intersection 
with the nodal lines of the d x 2_ y 2 gap function (see Fig. ||a). 

Under the above conditions one can easily rewrite the low-frequency asymptotics of 
damping by substituting ep — > t in the parenthesis and Bp — ► fi in other parts of Eq. (9). 
The range of validity for such expression is restricted by the condition uj <C min (A 2 /fi, A). 
Next, in order to obtain expression for damping at min (A 2 / /x, A) <C uj <C max (/x, A), the 
same substitutions in Eq. (10) should be done. Obviously, all the statements concerning 
the areas of principle contributions still hold. Finally, at max (/x, A) <C uj <C t we obtain 



Thus, the additional linear-in-a; high-frequency regime appears in this case. It occurs 
since at to ^> max (/x, A) the areas of integration in Eq. (5) engage the saddle-points of the 
normal-state spectrum located at (±7r, 0) and (0, ±7r). For the normal state crossover to the 
linear behavior has been demonstrated in ||. 

It follows from the aforementioned that there is no room for linear regime of damping 
at fi ~ t while the crossover from the cubic dependence to the quadratic one occurs at 
u ~ A 2 / fi. This case can be treated as the one of slightly deformed circular FS without 
additional characteristic scales. When fi decreases (i.e., close to the half-filling) the long 
almost flat pieces of the FS appear as well as the saddle-points in the normal-state spectrum 




(11) 
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come into play. For /i < A < t there co-exist three ranges with cubic, quadratic and 
linear frequency dependencies matching at uj ~ A 2 //i and u) ~ /x, respectively. At last, for 
very small /i<C A linear approximation for spectrum works well unless the ellipses touch the 
saddle-points. Consequently, the range of quadratic dependence in the latter case disappears 
and there exist the only crossover scale u ~ A. 



C. Tight-binding model, d xy -wave pairing 

Finally, let us consider the situation when the nodal points lie at the points of max- 
imum of FS curvature. To illustrate this case we write the normal-state spectrum in 
the tight-binding form and choose the gap function as <i-wave with xy-symmetry, [[| 
A q = A sinq x smq y (see Fig. |b). Then the straightforward calculations at u <C /x yield 
the formulas similar to those obtained in previous subsection: 

/T/\ 2 ,,3 a2 

ImE(k n ,o;)~-(-J — , "<-£r- (12) 



When A 2 n/t 2 C u C /i we have 

' uj t 2 fi * 



'V\ 2 uj 2 
ImEfkn,^ ~ - ( - —In 



t J fx 



mm 



A 2 /i ' uj 



(13) 



Again, we obtain the crossover from the cubic regime to the quadratic one. We see that 
close to the half-filling (at /i<f) the energy scale of this crossover (A 2 /i/t 2 ) is much less 
than that previously considered. It is connected with large curvature of FS near the nodal 
points so that the linear approximation for the spectrum breaks down much earlier. Note 
that this narrowing of the uj 3 region is accompanied by the strong parametric enhancement 
of ImS by factor t 2 /fi 2 . 

At higher frequencies /i <C w C t we get 

. f uj t 



V\ 2 

ImSfka, £j) ~ — ( — ) uj In 



mm — ; — 

\H uj j 



(14) 



Comparing this formula with formula found in previous subsection we see that the additional 
logarithmic factor exists stemming from closeness of saddle-point to the nodal one. 
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Note also that under /j->0we deal with poor-defined excitations since as it follows from 
our calculations ImS ~ ulnuj for all uj in this case. From Eq. (14) and Kramers-Kronig 
relations one can find that the quasiparticles are well-defined if 

2 



t fj, 



< 1 



(15) 



which is nothing but the criterion of absence of spin density wave. We note here that the 
same criterion arises in the case considered in Subsect. B, as well. The additional logarithmic 
enhancement stems in that case from the self-energy behavior near the "corners" of FS (cf. 

!)■ 

IV. DISCUSSION AND CONCLUSIONS 



All the results presented in previous section may be rewritten in an universal form. 
To do that more compactly let us introduce the following variables a = v 2 /(vplC) and 
b = min (v g , v g /JC), where /C is the curvature of FS at the nodal point. We define also a as 
the angle between vp and v 9 . Then we may write down the following formula: 



ImS(k„, uj) 



V 2 



uj <C min (a, b); 



v F Vg sin 2 a 



x < 



r . (u 


b 2 \ 


mm — , 




\ a 


auj J 



ui 2 a In 



ub 2 x ( 1 or In 



au) 



( b 2s 

min (a, b) uo max b, — 

— < OJ, 
a 



(16) 



This result is valid until v g <C vf and a ^ 0, as it takes place in all principal cases. The 
additional extrafactor In (aui/b 2 ) at high frequencies uj ^> b 2 /a appears only when the nodal 
point lies near the saddle one. It is easy to see that all previous results are reproduced by 
the following substitutions (i) /C ~ 1, v g ~ A, vp ~ Ef (assuming qp ~ 1) for circular 
FS, (ii) JC ~ fi/t, v g ~ A, vp ~ t and (iii) /C ~ \ftjjj>, v g ~ A^JI/t, vf ~ y/Jd for FS of 
varying curvature. We mention also that the quadratic regime disappears completely when 
/C v g /vF and that the intermediate crossover scales differ principally as K < 1 and K. > 1. 



Note that we have evaluated the damping using spectrum expansion near specific points 
up to the second order terms in deviations from this points. This is obviously valid for Fermi 
surfaces smooth enough. 

We undertake our consideration for circular and simplest tight-binding FSs. Meanwhile, 
in the adequate theory of high-temperature superconductivity the next-nearest-neighbor 
hopping terms in the tight-binding scheme should be taken into account. |10| It is clear from 



our analysis how the picture will be modified in this case. Say, for <i x 2_ y 2-wave pairing the 
range of cubic regime will be shrunken and a set of intermediate crossovers may appear in the 
range of the quadratic one. However, they will differ only by the logarithmic factor. Next, 
for the special choice of the parameters 2t' ~ t and /x ~ At' the saddle-points near (0, ±7r) 
and (±7r, 0) appear to be the "extended" ones (quasi-lD). [[11] Here t' is the next-to-nearest 
neighbor hopping integral and /i controls the filling. This case should be treated separately. 

The above consideration was made for zero temperature. It can be trivially generalized 
onto low finite temperatures (less than uj and minimal frequency scale in (16)). In the 
opposite limiting case (but not very close to T c ) Eqs. 8-14 can be rewritten by substitution 
T for u. In particular, it leads to the cubic T-dependence of damping for sufficiently small 
temperatures: u <^T <^ min (a, b) (the temperature dependence of the gap function yields 
only the corrections for this law). In order to apply our results to the normal state one should 
put simply v g = 0. Note that divergence in the low-temperature asymptotics never exist, 
since always T c 3> a. If the crossover scales min (a, b) or max (b, b 2 /a) are found below T c , T 2 
or T regimes should appear in the superconducting phase. When b 2 /a ^> T c , the crossover 
from the conventional 2D FL behavior to the linear one occurs in the normal state. 

Our treatment has demonstrated, in particular, that the "purely nodal" regime for damp- 
ing takes place in the quite narrow frequency (temperature) range. It is restricted from above 
by the much less scale than the gap value as one could expect. It seems to be not a feature of 
damping only, so that the same characteristic scale should appear for other quantities such 
as density of states, nuclear relaxation rate, penetration depth, inelastic scattered neutron 
intensity near incommensurate and internodal vectors, etc. Since, on the other hand, the 
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impurities scattering restricts this range from below, one can estimate it as follows: 



where T imp is the damping due to impurities. The estimates show that the range defined by 
Eq. (17) appear to be quite narrow for the high-T c compounds. It may even disappear if 
Timp will be large enough (~10mel / ). 

In conclusion, we consider the damping of excitations in the 2D superconductor with 
nodal points in the spectrum. We show that at T = this damping demonstrates gap- 
less cubic-in-cu low-frequency behavior at the nodal points. The results obtained for high 
frequencies depend on details of Fermi surface parametrization. We have analyzed two im- 
portant particular FS parametrizations such as the models of almost free and tight-binding 
electrons. We find the frequency variation of damping and the crossover scales for these 
cases. For smooth enough FS they may be expressed through the model-free characteris- 
tics of the spectrum at the nodal points. We argue that the excitations are well-defined in 
the wide frequency range. We extend our treatment to finite temperatures and study the 
possibility of crossovers in temperature dependence of damping below the superconducting 
transition. 
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APPENDIX A: 



Here we study Eq. (5) for circular FS: £ q = q 2 /2m — ep. We shall assume that Vk,q; P ,i is 
finite in the important areas of the integrand in Eq. (5). Also, let us take for a moment the 
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combination of u — v coefficients in Eq. (5) to be of order of unity. Then we should analyze 
the following expression 

ImE ~ V 2 jj D d 2 qd 2 p 5 ( £q + e p - w + £ q + P - k ), (Al) 

where D is defined by the condition 

£ q + e P < u- (A2) 

It is convenient to make the substitution q = q 2 — q 2 F , p = p 2 — q 2 F in Eq. (Al), then 
d 2 q = dipdq, d 2 p = dipdp, where ip and ip are the corresponding polar angles. Taking q F ~ 1, 
we may represent the normal-state dispersion as £ q ~ e F q. Besides, we expand the gap 
functions up to the linear in angle terms A q ~ Aip, A p ~ Aip. Using such parametrization 
and leaving the principal terms only, we reduce Eq. (Al) to the following one 

min(l,w/A) u/ep 

ImE -V 2 J j d(pdi> J J dqdpb (V<? 2 4 + A V 
o o 

+s]p 2 e 2 F + A 2 ^ 2 -u + e F <pi?j , (A3) 

where the limits of integration stem from inequality (A2). This condition means, in par- 
ticular, that at small u momenta integrations in Eq. (A3) are restricted by the elliptic 
areas whose centers lie at the nodal points. It is convenient to pass therefore to the elliptic 
coordinates in Eq. (A3): 

r = yV4 + A V R = \jv 2 z\ + A V 
a = arccosge^/r (3 = &rccospe F /R 



(A4) 



Then it transforms to 



min(l,A/w) 

V 2 u 6 



ImE ~ 2 ^ 2 J J da dp J J dr dR 



■ F o 



.'r + R—le F 
xS[ + — ua(3) (A5) 



rR A 2 

This expression allows us to analyze all principal limiting cases. 

(1) uo <C A 2 /e F . In this case we can neglect the second term in 5-function. The integral 
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remained can be easily evaluated and the result is given by Eq. (9) 

(2) A 2 /e F CwCA. Here it is convenient to rewrite 5-function in Eq. (A5) as follows 



One can see that the integration over one of the angles gives the logarithmic divergence 
which is cut on the scale A 2 /ueF- 

(3) A <C uo <C e F . Then we should integrate the 5-function in Eq. (A5) with respect to the 
angles between the limits and A/a;. After the trivial substitution (a, (3) h- > (A/a;) (a, (3), 
we have 



Again, we have the logarithmic divergence, which is cut now on the scale uj/sf- The last 
two cases we may be combined as Eq. (10). 

It can be checked out that taking the combination of u — v coefficients in Eq. (5) into 
account does not change our estimates. 

Other kinds of Fermi surface and gap functions considered in the text can be analyzed 
on the similar manner. 




(A6) 




(A7) 
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FIGURES 

FIG. 1. The disposition of the momenta which give the principal contribution to ImS at u — > 0. 
The solid lines denote the Fermi surface and the dashed lines correspond to the nodes of the gap 
function. 

FIG. 2. The disposition of the momenta which give the principal contribution to ImS at to 3> \i 
in the case of nested FS. The crosses label the saddle-points. The nodal points are located (a) far 
away and (b) near the saddle ones. 
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